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1 Problem 

Magnetostatics can be defined as the regime in which the magnetic fields B and H have 
no time dependence, and "of course" the electric fields D and E have no time dependence 
either. In this case, the divergence of the fourth Maxwell equation, 

4:71 1 <9D 

V x H= — J free + -— , 1 
c c at 

(in Gaussian units) implies that 

V ■ J frcc = 0, (2) 

i.e., that the free currents flow in closed loops. Likewise, the time derivative of the fourth 
Maxwell equation implies that Jf ree has no time dependence in magnetostatics. 

Often, magnetostatics is taken to be the situation in which V • Jf rec = and D, E and 
Jf ree have no time dependence, without explicit assumption that B and H also have no time 
dependence. Discuss the possibility of waves of B and H, consistent with the latter definition 
of magnetostatics [T]. 

Consider two specific examples of "magnetostatic" waves in which Jf rce = 0: 

1. Ferromagnetic spin waves in a medium subject to zero external field, but which has 
a uniform static magnetization that is large compared to that of the wave. That is, 
M = M z + me ! ' k ' r ~ ui ', where m <C M . Here, the quantum mechanical exchange 
interaction is the dominant self interaction of the wave, which leads to an effective 
magnetic field in the sample given by B c ff = aV 2 m, where a is a constant of the 
medium. 

2. Waves in a ferrite cylinder in a uniform external magnetic field parallel to its axis, 
supposing the spatial variation of the wave is slight, so the exchange interaction may 
be ignored. Again, the time-depedendent part of the magnetization is assumed small 
compared to the static part. Show that the waves consist of transverse, magnetostatic 
fields that rotate with a "resonant" angular velocity about the axis. 

In practice, the spin waves are usually excited by an external rf field, which is to be neglected 
here. 

2 Solution 

2.1 General Remarks 

In both definitions of magnetostatics the electric field E has no time dependence, d~Ei/dt = 0, 
so the magnetic field B obeys d 2 B/dt 2 = 0, as follows on taking the time derivative of 
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Faraday's law, 

(in Gaussian units). In principle, this is consistent with a magnetic field that varies linearly 
with time, B(r, t) = Bo(r) + Bi(r)£. However, this leads to arbitrarily large magnetic fields 
at early and late times, and is excluded on physical grounds. Hence, any magnetic field B 
that coexists with only static electric fields is also static. 

There remains the possibility of a "magnetostatic wave" in a magnetic medium that 
involves the magnetic field H wavc and magnetization density M wavG which are related by 

= B = H + AirM. (A) 

w - LJ wave -"--"-wave ~ ^* " - LVJ -wave • V / 

If there are no free currents in the medium, and any electric field is static, then the fourth 
Maxwell equation is simply 

V x H = 0, (5) 
which defines a subset of magnetostatic phenomena. 



2.2 Ferromagnetic Spin Waves 

Consider a ferromagnetic material that consists of a single macroscopic domain with mag- 
netization density M = Mqz + m(r, t), where Mq is constant and m <C Mq. We suppose 
there are no external electromagnetic fields. Associated with the magnetization M are mag- 
netic fields B and H whose values depend on the geometry of the sample. We suppose 
that the weak time-dependent magnetic fields due to m lead to even weaker time-dependent 
electric fields, such that the situation is essentially magnetostatic. The consistency of this 
assumption will be confirmed at the end of the analysis. 

The ferromagnetism is due to electron spins, whose dominant interaction is the quantum 
mechanical exchange interaction, in the absence of external fields. For a weak perturbation m 
of the magnetization, the exchange interaction preserves the magnitude of the magnetization, 
so its time evolution has the form of a precession j2] , 

^ = dxM. (6) 

dt v ; 

As this is the same form as the precession of a magnetic moment in an external magnetic 
field the precession vector Q is often written as a gyromagnetic factor T = e/2m e c ~ 10 7 
Hz/gauss times an effective magnetic field B e g (or H e ff). Here, e > and m e are the charge 
and mass of the electron, and c is the speed of light. For a weak perturbation in an isotropic 
medium 

B cfr = «V 2 m, (7) 

where a is a constant of the medium. 

Then, the equation of motion of the magnetization m is 



dm 

~dt 



oTV 2 m x M. 
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For a plane- wave perturbation, whose phase factor is e l ^' r ~ ut \ the equation of motion 
(jBJ) becomes 

icum = aTk 2 m x M z. (9) 
This is satisfied by a circularly polarized wave, 

m = m(x + iy)e i(k - r - w ' ) , (10) 

that obeys the quadratic dispersion relation jl] 

to = aTM k 2 , (11) 

which implies that u <C ck in physical materials, where c is the speed of light. Hence, the 
electric fields are much smaller than the magnetic fields associated with the time-dependent 
magnetization m, so that V x H = to a good approximation, and we may use the term 
"magnetostatic" to describe the waves. These waves of magnetization are, however, better 
termed "spin waves" , whose quanta are called "magnons" . 



2.3 Rotating Magnetostatic Modes in a Ferrite Cylinder 

In the magnetostatic approximation the fields B and H obey 

V-B = V-(H + 4vrM) = 0, V x H = 0, (12) 

where the field B but not H and M must be static (or at least so slowly varying in time that 
the resulting electric field is small compared to B). We first consider a ferrite of arbitrary 
shape of characteristic length a in a uniform external magnetic field B cxt = H cxt = Hqz. We 
suppose that this field is strong enough to induce a uniform magnetization Mqz throughout 
the sample. 

For waves with weak spatial dependence as we shall assume, the exchange interaction 
is negligible, since it varies as the second spatial derivative of M. Then, the spins interact 
primarily with the local magnetic field B according to 

^ = rBxM = rHxM. (13) 
dt v ; 

We consider a perturbation m to the magnetization that has frequency uj and wavelength 
large compared to the size of the the sample. Then the total magnetization can be written 

M = M z + me"", (14) 

where m -C M . Similarly, we write the magnetic field inside the sample as 

B = B z z + be~ iwt , H = H z z + he~ iuJ \ (15) 

where B z = H z + 4irM and H z = H Q — 4ttN z M are the sum of the external field and 
that due to the uniform magnetization Moz, and so are also uniform for spheroidal (and 
cylindrical) samples whose axis is the z axis [3]. The "demagnetization" factor N z varies 
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between 1 for a disk and for a cylinder. The perturbation m exists only inside the sample, 
but the corresponding perturbations b and h exist outside the sample as well. 

Inserting eqs. (|14|) and (|15|) in the equation of motion (|13|). we keep only the first-order 
terms to find 

-ium = Fz x (M h- H z m), (16) 



whose components are 



r 

m x = i—(M h y -H z m y ), 

r 

m y = -i—(M h x -H z m x ), (17) 
m z = 0. 



(18) 



We solve for m in terms of h as 



m x = ah x - i(3h yj 

m y = i(3h x + ah y , (19) 

where 

_ T 2 H Z M _ TM uj 

For later use, we note that in cylindrical coordinates, (r, 9,z), eq. ()19|) becomes 

m r = ah r — i(3hg, 

rag = i(3h r + ahg. (21) 
As we are working in the magnetostatic limit (|T2*|) . we also have 

V ■ b = V ■ (h + 47rm) = 0, V x h = 0. (22) 
Hence, the perturbation h can be derived from a scalar potential, 

h = -V0, (23) 

and so, 

V 2 = 4ttV • m. (24) 
Outside the sample the potential obeys Laplace's equation, 

V 2 = O (outside), (25) 

while inside the sample we find, using eq. (|19|). 

(1 + 47ra) ( -d + w ) + d$ = (inside) ' (26) 
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The case of an oblate or prolate spheroid with axis along the external field has been 
solved with great virtuosity by Walker jSj, following the realization that higher-order modes 
deserved discussion jjj. Here, we content ourselves with the much simpler case of a long 
cylinder whose axis is along the external field, for which the lowest-order spatial mode was 
first discussed by Kittel jHj. We consider only the case of waves with no spatial dependence 
along the axis of the cylinder. 

With these restrictions, both eqs. (J25j) and (J26)) reduce to Laplace's equation in two 
dimensions. We can now work in a cylindrical coordinate system (r, 9, z), where appropriate 
2-D solutions to Laplace's equation have the form 



(r<a,9) = J2 r -(A n e me + B n e- me ), (27) 
or 

n 

(r>a,9) = J2^(A n e me + B n e~ inff ), (28) 



which is finite at r = and oo, has period 2ir in 9, and is continuous at the boundary r = a. 

The boundary conditions at r = a in the magnetostatic limit (|22[) are that b r and hg are 
continuous. The latter condition is already satisfied, since hg = —{\/r)d(j)/d9. We note that 

b r = h r + 47rm r = (1 + 47ra)h r - Ani(3h e , (29) 

recalling eq. 1)21)1 . Using eqs. (J27)) and (J28)) we find that continuity of b r at r = a requires 

77 r 

Y - (1 + 2vra + 2nP)A n e me + (1 + 2vra - 2Ttp)B n e~ me = 0. (30) 



Nontrivial solutions are possible only if 2ir(a ± (3) = —1, in either of which case there is an 
infinite set of modes that are degenerate in frequency. Using eq. ()2"U)) . we find the "resonance" 
frequency to be 

u = ±T(H + 27iM ), (31) 

noting that for a cylinder the demagnetization factor is N z = 0, so that H z = H , as is 
readily deduced by elementary arguments. Since we consider frequency u to be positive, we 
see that the two solutions ()3*T)) correspond to two signs of Hq, and are essentially identical. 

For spheroidal samples, the modes are enumerated with two integer indices, and are not 
all degenerate in frequency, as discussed in j^J. 

We close our discussion by showing that the electric field of the wave is much smaller 
than the magnetic field. The scalar potential for mode n is 

Mr <a) = r ^e l{n6 -^\ K{r > a) = ^"). (32) 

We see that for n > the potential rotates with angular velocity Q n = u/n about the z axis. 
The potential is maximal at r = a, so consistency with special relativity requires that 

auj 

v(r = a) = — c, (33) 
n 

which appears to have been (barely) satisfied in typical experiments jH]. We also see that 
for high mode number the spatial variation of the wave becomes rapid, and the neglect of 
the exchange interaction is no longer justified. 
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The magnetic field h = — V0 of mode n has components 

hr(r <a) = - n r -^-e^ 6 -^ = ~(j> n , h r (r > a) = n^e^ 9 '^ = (34) 

hg(r < a) = ih r (r < a), hg(r > a) = —ih r (r > a). (35) 

The monopole mode, n = 0, does not exist. The lowest mode is n = 1, which corresponds 
to a uniform, transverse field h that rotates about the z axis with angular velocity u>. 
From eq. (|21|) we find the magnetization to be 

m — (36) 

Z7T 

for all modes (for r < a only, of course), so the magnetization of mode n also rotates with 
angular velocity u/n. 

The magnetic field b = h + 47rm is then, 

b(r < a) = — h(r < a), b(r > a) = h(r > a). (37) 

Using either the r or 6 component of Faraday's law, we find that the associated electric field 
e has only a z component, 

e=^0z, (38) 

both inside and outside the cylinder (consistent with continuity of the tangential component 
of the electric field at a boundary, and with V • e = 0). The ratio of the electric to the 
magnetic field of mode n at r = a is 



— , (39) 
nc 



which is small so long as condition (|33J) is satisfied. Hence, the condition that auo/c 1 is 
doubly necessary for the validity of this analysis. 

Because the magnetization m is moving (rotating), there is an associated electric polar- 
ization p according to special relativity 



v 

7— x m. (40) 

c 



For mode n we have v = uirO/n <c, so7 = l/yl — v 2 /c 2 ~ 1, and 

urm r A urh r A cu . A e 

p = z = z = <p z = . (41) 

nc 2irnc 2nc 2n 

The electric displacement is related by d = e + 47rp, which has the value 

d(r < a) = — e = ——0 z , d(r > a) = e = —0 z , (42) 
c c 

The fourth Maxwell equation now implies that 

1 <9d lo 2 

Vxh=-Vx V0 = -— = ±i— 0z (43) 

C OT & 

Thus, there is a small violation of the magnetostatic conditions (j22j) . but this is second order 
in the small quantity au/c (noting that V0 ~ <f>/a). 
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